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Let  Xj - XN  be  a  sample  of  N  Independent  observations  with  a 

normal  distribution  N(m,cr-).  In  1935,  W.  R.  Thompson  has  obtained 
the  probability  distribution  of  an  observation  X?  chosen  at  random 
in  the  sample,  when  the  mean  m  and  the  standard  deviation  cr  are 
replaced  by  the  correspond  I ng  sample  character? st ics  [I];  namely 
(Xj-302/{N-J)S2  follows  an  Incomplete  beta  distribution  with 
parameters  1/2  and  (N/2)-!. 

This  result  has  been  generalized  in  several  directions,  in 
[2];  the  distribution  of  a  sample  is  given  after  its  being 
centered  and  student ized  by  means  of  the  mean  and  standard  devia¬ 
tion  of  an  independent  sample;  the  distribution  of  a  subsample 
is  derived  when  m  and  or'  are  replaced  by  the  corresponding  sample 
character  1st ics,  and  the  conditional  distribution  of  a  subsample 
is  obtained  given  7  and  S. 

Further  generalizations  have  been  presented  in  [3]  where  the 
results  just  mentioned  are  extended  to  the  case  of  the  multivariate 
normal  distribution.  Applications  to  Bombing  Theory  were  proposed 
in  [4] o 


2 


I.  I 


Let  (^  )  =  £. *  with  =  (X  j , .  .  ,X^, .  .Xf)  be 

a  subset  of  observations  from  a  set  of  M  Independent  observations 
(X)  <=  (X|,.0,Xj,« . X^j) *  «*  (Xp  with  the  same  p-variate  normal  dis¬ 
tribution  N(m,^);  let  m*  and  S  be  the  vector  mean  and  covariance 
matrix  of  the  subset  The  distribution  ofSqfven  m*  and  S 

for  k  4  M-p  is 


?(T5|m»S)d?  -  C.|l  -  (k/N)S-,S<  -  [k/(H-k)3S-'S*| 
in  the  domain  where  the  determinants  ore  positive,  with 
S'"  =  (m^  -  m’ )  (irl^  -  in’) 
and 

-  r'r'  r  *'  "'■»  r'rty  ,  P[(N-k-i)/2j..,Pl C(M-k-p;/2j. . * 

^kp/2N(l;-l)p/2(|]_k)p/2 
This  result  was  given  in  C 3 U *• 

In  view  of  the  completeness  of  m’S,  it  results  from  Blackwell’s 
theorem  that 

j...  [if  (5  )  f  |  m’S)d? 

will  provide  a  uniformly  minimum  variance  unbiased  estimate  of 
G[  (  5  )  ]  if  it  exists. 

As  a  special  application,  if 
g(A,m,^t)  =  j  ...  |lUm,^)dX, 

a  minimum  variance  unbiased  estimate  of  g(A,m,  i>  will  be  g i ven  by 

\  ...  \  i  (5  j  m*  S)d  £  with  k  =  !, 

U  > 

since  an  unbiased  estimate  ot  g(A,m,>_ )  is  provided  by  the  charac¬ 
teristic  function  I  (?  )  of  the  set  A  whose  conditional  expected 

A 

value  is  given  by  the  integral  above.  Further,  in  view  of  the 
fact  that  summation  keeps  unbiasedness,  the  minimum  variance  un~ 


1 1 


(N-k-p~2)/2 


d?/  S 


k/2 


2.1 


biased  estimate  of  the  density  N(rn,  ^.)  is  f  (S 


m’S) 


-  3  - 


The  already  established  results  for  the  univariate  distribution 
can  be  obtained  by  considering  if  as  a  special  case  of  the  multi- 
var iate  si tuaf ion. 

2.2  The  conditional  probability  distribution  of  ^  given  m',S 

is  obtained  by  writing  the  joint  distribution  of  tne  mean  m'  and  co- 

n 

variance  matrix  S  of  a  sample  of  size  n  =  N-k  and  of  ml  and  S  , 

n  5  * 

then  performing  the  change  of  variable 


m*  =  Nm’/n  -  km’  /n 
n  5 

S  =  MS/n  -  kS.  /n  -  S*kll/n2 
n  s 

subsequently  dividing  by  the  distribution  of  S. 

m’  is  M(m,!j£/n)  distributed,  5^  is  '7(^1,  n- I )  distributed,  the 
Jacobian  of  the  change  of  variables  is 

DK»S">  |M,  ,p(p+3)/2 

-  _  (f]yn) 


D(m! ,3) 
m’  is  ll(m,^/N(, 


S  is  N-i)  distributed. 

After  grouping  terms,  oneob tains  for  p^k^N-p 


ft 


nlj,  \ 


m* , S)dm’ dS  = 

,  I  ( N-k-p-2 ) /2 

I-(k/N)S,  S“5  ~[k/(N-k)  jS-'S”  1 
1/2  ,  ,(p+l)/2  ' 


s5  s 


(k-p-2)/2 


( dm’  / 


S  )  dS5  /  5j 


where 
C  = 


L( N-k) N 


k-  I 


p/2  p(p+!)/4  p 

[l/TT  3TT  P  [N-j}/2] 

I 


2.3 


rcN.k.j)/2irco«-j)/23 

in  the  domain  where  the  determinants  are  positive., 

In  view  of  the  completeness  of  m*,S.  the  kernel 


f  (m 


t  s 


m*.S)  is  of  great  interest 


If  Sn)  is  an  unbiased  estimate  of  a  function  g(m,^£) 

of  the  parameters,  then,  as  a  consequence  of  the  uniqueness  of 


~  h  - 


minimum  variance  unbiased  linear  estimates  in  case  of  completeness 

■[ . 1  4Vm|'S5  )f<mPst 


rn*  ,S)dm*  dS? 


i,e,  minimum  variance  unbiased  linear  estimates  are  a  subclass  of 
the  solutions  of  this  integral  equation  which  defines  a  multivariate 
fransf ormat ion . 

3,0  In  case  X  is  normally  distributed  NCC^cr),  univariate  situa¬ 

tion,  or  N(0,1^),  multivariate  situation,  other  generalizations  are 
of  ? n t er es t . 

3,5  In  the  univariate  case,  let  X  be  a  set  of  N  observations 

S  ...Xjj  with  probability  distribution  N(0,  er ) ,  the  conditional 
distribution  of  a  subsarnple  ^  of  k  items  drawn  from  X  at  random, 
given  the  moment  of  order  two  rn^(X)  of  the  set  X,  is  obtained  by 
writing  the  joint  distribution  of  the  moment  of  order  two,  m^(Xn) 
of  a  sample  Xn  of  size  n  =  N-k  and  of  an  independent  sample  \  of 
k  items,  performing  the  change  of  variables 

rrd  ( X n )  =  (M/n)m*(X)  -  (k/n)m£(?  )  whose  Jacobian  is  N/n,  then 
dividing  by  the  distribution  of  m£(X),  one  obtains  with  ^ 


(X  ...Xk),  k 4  M 
ft?  |  m£(X)]dl  = 

P  ( N/2 ) 


k/2 

TT  p[(  N-k)/23 


l-Ck/M)  m*(X) 
2 


(N-k-2)/2 


dS 


CNm|(X)j 


k/2 


in  the  domain  where  the  bracket  is  positive. 


if 


•*-«//  ml(X) ,  one  obtains 


*m 


. (,  .xv  )(H-k-2)/2(dT/„k/z) 

V  rC(M-k)/2]  M  " 


'X’K  Z  H 


f.e.  the  distribution  of  a  subsampie  normalized  byym*(X)  is  de- 
pendent  from  m£. 

Also,  in  case  k  =  N- 1 ,  u  -  ?//  Nm|(X),  u j . . . ^  are  N-l  In¬ 
dependent  coordinates  of  a  N  dimensional  unit  vector,  and  one  has 
n  M/2  ,  -1/2 

Cl  (N/2)/ir  3(1  -  uu  )  du  uu  /.  I. 

The  distribution  of  any  coordinate  uj  being 


- ,  ( I  — u  ) 

/N~  I  l  v  { 


?  (M-3)/2 


3  i  ) 


2 

u  L.  I 

l  ~ 


(if  M  =  3  one  obtains  a  rectangular  distribution). 

Clearly,  the  quantity  v  ~\}  N-l  u.  i\/  l-u?  is  Student  distributed 
wit,:  N—  I  degrees  of  freedom. 

Uj  is  the  cosine  of  the  angle  4*  j  °f  ihe  unit  vector  with  any 
direction,  uj  =  cos  ip  •  and  v  =V n-l  cot.  Lp  j . 

This  result  can  be  slightly  generalized  to  the  case  of  a 
spherical  di str ? but  ion. 

Let  X  =  (Xj,...,Xn)  have  probability  distribution  f (X)dX, 
letR2=XX’,  let  if  (X)dX  be  spherical,  i.e.  f(X)*h(XX»)» 
then  R2  is  a  sufficient  statistic  for  f (X)dX  . 

f  (X)dX  =  h(R2)dX  =  h(R2)Rn_'ldRdXn  ,  (Scathe  area  of  the  unit 
sphere  in  the  n  dimensional  space);  i  nt  egra  t  i  ng  out  d2»n  j,  one 
obtains  5Zn,h(R2)Rn-^dR  as  the  probability  distribution  of  R,  a 
well  known  result. 


More  explicitly!  using  polar  coordinates 


X  j  =  Rcos  tp  j 

Xn_j  =  Rsinifj  sin  ipg  . . .sin ip  n_2costf  n_, 

=  Rsinif  j-,...  .esin  ip  n^2s  in  tp  n_j 
Vp  .  in  ( o , ir  )  if  i  4  n- 1 ,  vf  j  i n ( o ; 2 tt  )  if  i  =  n~ ! 
dX  =  JdRd  Ip  |o.dlfn_3,  where 

J  =  [D(X)/D(R3vp  )]  =  Rn"!sinn  2cpj  .  .  .si  n  ip  n<_2 

and  dZ.  nl  =  sinn~2  ip  j  » .  .sin  tp  n_2d  if  p  «  »d  }  =  JdvP  /rD  * 

=  dX,  ...c3Xn_]/|xn|Rn-2 

hence;the  conditional  distribution  of  X  given  R  is 


d  21  IT|.  —  P  (  n/2) 

~~Z  ™  n?2  s,n  po.sinlf  n„2d  H*  |°',adl?  n„s 

t-  n  5  2  7 T 


it  is  the  distribution  of  the  independent  polar  coordinates  of  a 
unit  vector  equ id istr i buted  on  the  unit  sphere.  In  terms  of  the 
Xj’s  and  u.’s,  this  gives  (taking  into  account  both  signs  of  Xn), 

P  ( n/2)  dX,,..dXn_i  -i/2 
-^75-  =  du 


Consequently  the  distribution  of  any  subsample  is  the  one 
given  in  paragraph  3-1. 

This  shows  that  "Thompson’s  distribution  depends  only  on  the 
sphericity  of  the  universe.  If  we  consider  a  random  unit  vector, 
with  fixed  origin,  "chosen  at  random"  in  the  n  space  i.e.,  with  a 
probability  d istr i but  ion  invariant  under  the  rotation  group,  we 
.  see  that  we  can  construct  such  a  vector  by  normalizing  a  vector 
whose  distribution  is  spherical. 


Let  us  consider  a  set  (X)  of  independent  N  observations 


of  o  p-rd  imensi  ono  i  vector  with  distribution  N(oa  ?  )»  (X)  = 

=  (xf)  =  (X’s...,xJ  „.XP)  XJ  =  (XJ*»eoX^0„XJ)' 

*  t  ‘  '  j  N 

let  (i  )  be  a  subset  of  K  observations-  Let  S  and  be  the 

maximum  likelihood  estimates  of  ^  obtained  respectively  with 

(X)  and  (?).  The  distribution  of  (5)  given  the  sufficient 

statistic  S,  is  obtained  by  writing  the  joint  distribution  of 


{?  )  and  the  maximum  likelihood  estimate  S  of  %  obtained  with 

n 

a  sample  of  n  ~  M-k  observations,  making  the  change  of  variables 
3n  *  (M/n)  S  -  (k./n)S1 

end  dividino  by  the  distribution  of  S*One  obtains 
/>.  I  (N-k-p- 1  )/2  ,  ,  k/2 

Jf(5|s)d?  =  Co  j  :c  -  (k/ins  s%  d?/  |s 


c  -  T7-pC(H-H-j).- 13  pk  •' 2 

V _ _  [  •  /  (  N  7r  )  j 

P'~ 

yj7‘»C  ( N-k-H  -  j)  i  2] 

vs)  id  for  psN-k,  in  the  domain  where  the  determinants  are  positive 
Now  kS,_  =  ?  £  j 

so  •;  he  r  k  k 

I  -  (k/N)S‘  S?  |  =  I  -(s“i/H)Zl  S\  5  =  |ns|'-,|  S-EIj  t 

there  exists  a  unique  triangular  matrix  T  with  positive  diagonal 
such  That  NS  -  X' X  =  T  T*  then 

I  -  (k/M)S*fS.  I  =  [  T  S-ZS'£.  T' 

s  !  >  j  J 

i-  2  t-'s’It'-'  -Ii  -  Z  V*J  if 


8 


Mow,  it  is  known  that 

i-^pj  -Ki t-Vrft’ 


so  that 


—  1, 


I  -  (k/N)S~ 

-f  (  F  S)d  £  =  C.J  t)~  -  X  j  (S-’/N)^ 


K  -  ^iT’-'-r'St  |  -  |^.f  -  5,<s-' 


if  >! 
(N-k-p-!)/2 


dF  /  S 


/  N )  ^ 

k/2 


One  can  a  iso  work  with  the  variables  making  the  transforrnat  ion 
£  T*  ”  *  ,  one  obtains  for  the  distribution  of  ■??  , 

(N-k-p-l )/2 

d 

k  i  |  k/2 


c: 


1  -  9  'y 


since  tne 


^  .  |k  |  k/2 

Jacobian  D(  F  )/D{'y  )  =  I  TJ  =  SJ  , 

7“  —  j 

*)  =  feT* 

is  a  procedure  of  si  ud  enf  5  za  t  i  on  of  the  set  ?  and  generalizes  the 
student  izat  ion  by  R  =y  Nm^(X )  in  the  univariate  case  to  the  multi- 
variate  case. 

in  some  problems  the  fol  towing  remark  is  of  interest* 

Consider  theor  thogona  I  matrix  0  that  diagonalizes  (S"vN)  into 

3  t 

f\~  say,  then  pr  emu  I  f  i  p  I  y  i  ng  and  pos  tmu  1 1  i  p  I  y  i  ng  by  0”'  and  0  we 


have 


U  '  (  k/N)S~*S 


A 


- 1 


o“!  F »  f;  .o 

j'  j 


=  |  r  -  A 


~  i. 


u!u 


Pi 


where  u.  is  the  vector  of  the  coordinate  of  F  .  vv.r.  to  the  system 

of  coordinates  constituted  by  the  eigen  vectors  of  (S"*/N)  i.e,  of 

NS,  (which  is  almost  surely  of  rank  p). 

One  can  consider  the  norma  1  i  zat  i on  of  these  coordinates  by 
. _  |  /2 

the  roots  yhj  •  •  yhp  of  f\  3  we  get,  in  terms  of  the  new  vari¬ 


ables  V 


:»• 


-  9  - 


,  .  j ( N-k— p- i ) /2 

C  l  -  V  V  dV 


and  the  correspond ing  norma  I izat ion  is 

v-  Ss"-|/2 

|  /2 

according  to  the  usual  definition  of  S  (see,  for  Instance,  C  1 3U )  = 

The  conditional  distribution  of  given  S  Is  obtained 

by  writing  the  joint  distribution  of  and  Sn,  making  the  change 

of  variables  indicated  in  3»^  and  dividing  by  the  distribution  of 
S,  one  obtains 

D  |  I  ,  I  (N-k-P-l)/2  I  -,l (k-p-l )/2 

f  (S?  j  S)dS  =  C.  J  X  -  (k/N)S?  S  |  Is?  S  ’ 

,  |(p+})/2 


?T[(N+i~j)/2] _ 

jjp  C  ( k+i- j)  /2]  ft  ft  ( N-k+l- j)/2] 


dS?  /  js 
with 


C  =  (k/N)pk/2 _ I 


p  L  k  l.  N-p,  in  the  domain  where  the  determinants  are  posit ive0 
In  terms  of  4?  one  obtains 

*  j  (N-k~p-!)/2|  ,  j  (k-p-l)/2  , 

c  i-  *>">1  h  i 

» 

since  the  Jacobian  of  the  transformat  ions  from  ?  %  -  kS^  to 

y  is  T  =  J  S  J  .  This  is  a  generalization  of 

the  incomplete  beta  distribution. 

The  results  in  3.4  and  3-5  can  be  generalized  to  the  case 
of  a  probability  distribution  which  depends  on  the  set'  X  through 
X^X  =  NS.  ?„e„ 


f (X)dX  =  h(X*X)dX 


/ 


-  i 


-  JO  - 

Then  X’X  is  a  sufficient  statistic  for  the  distribution.  it  has 
been  shown  by  Hsu  that  the  probability  distribution  of  X*X  =  NS, 
under  such  circumstances,  is 

g (X’X) dX’X  =  7TNp/2  -P(P~!  )/**  |  x,x  ^ H  P  !  ^/2h(X*X)  dX*X  C93 

Tl.  fC(N-j+i)/2] 

I 

it  is  well  known  that  there  exists  a  unique  triangular  matrix  T 
with  positive  diagonal  such  that  X’X  =  TT*  andthat  the  transfor¬ 
mation  X  =  YT’  defines  the  socafied  "rectangu  I  ar  coordinates"  T 
[8]0  To  a  considerable  extent,  what  follows  overlaps  the  theory  of 
rectangular  coordinates,  though  the  approach  is  somewhat  different. 
3.6.1  The  vectors  (X1 , .  , .  ,X j, . .  „XP, ) ,  span  a  space  that  is  almost 
surely  p-d  rmensiona I .  'Ve  want  to  construct  an  orthonorma!  basis 
( Y*, . .  ,Y^. .  .YP)  in  that  space.  This  can  be  done  by  use  of  the 
Schmidt’s  orthogona fi zat i on  process,  starting  with  X*. 

-  U(yJ...,Y^...YP) 

where  U  is  an  upper  triangular  operator  (i.e.  in  the  system 
(Y|,..,YP),  u  is  represented  by  an  upper  firrangular  matrix  with 
posi t i ve  d iagona I . ) 

now  X’X  -  NS  =  ((xJxJ))  =  ((UYjUYJ))  =  ((yJu’UYJ)) 
i.e.  the  i,j  element  of  U’U  matrix  of  U’lJ  with  respect  to 
(yJ...,Yp)  is  the  i ,  j  element  of  NS  that  is  U’U  =  TT*  and 
from  the  uniqueness  of  the  factorization  U  =  T*;  consequently 
X  =  YT’ 

Y  has  orthoqonai  columns  and  Y’Y  =  J. 

,  ~  PP 

X  .  =  t  ;  Y'  +...+t j SYJ 
J  JJ 


X  involves  np  free  random  vari  abl  es,  Y,  np~p(p+I)/2  as  orthogonal 
and  T,  *  p(p+j)/2  as  triangular, 

i  _  | 

Y  =s  x  T  is  a  random  orthonormal  basts. 

Now 

»  • 

-P (X)dX  =  h ( X  *X )  TT  dXJ  „.dXJ 

j  N 

let  us  study 

TT  dXJ  ...dXj  »  TT  dXj 
j  1  N  j 

One  has  X^  =  (f.  Y*  +  ...  +  f +  fstY^ 
j..  J  J—  I  * 

Consider  the  j-  j  dimensional  space  spanned  by  Y^aCYd""l  and  complete 
that  basis  by  an  orthonorma!  basis  in  the  N-j+l  dimensional  comple¬ 
mentary  space;  i  . -YJ  is  the  projection  of  XJ  in  that  space  where 
J  J  i*  j* 

its  coordinates  are  XJ  jo.,,X  say.  in  the  new  N  dimensional 

J  N  j* 

basis  X~*  has  coordinaies  *'js!*o0^j  anc*  s‘nce  one  passes 

from  the  old  basis  to  the  new  basis  by  a  rotation  one  has 

dxd  =  dt;  ....df.  .  ,dXJ.  . .  0dX.d* 
d  >  1  J ,  J*- 1  J  » 


Using  polar  coordinates  in  the  N-j+J  complementary  space  we  have 

dX-l*.  ,,dXj*  =  fN”j  dt--dX  tl 

J  N  jJ  u  TJJU  ^l^j+S 

where  d£  .  denotes  the  elementary  area  of  the  unit  sphere  in 
N-J+J 

the  N-j+3  dimensional  space. 

Therefore, 
dX 


dT’ 


J  JJ 


since  dT'  =  2~p  TT1  j§+d  !dXX  and  Jf  t 
J  j 

(N-p-l )/2 


and 

dX  =»  2~P  I  X  *X 


JJ 


^dIN-j+|dX’X 


T 


one  has 


k 
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dF,  symbolizes  a  differential  expression  involving  N-j 
N-j+l 

angles  functionally  independent  [i.e.  Np-p(p+!)  variables  for  all 
j  from  I  to  p„] 

(N-p-0/2 


i  (X)dX  =  h(X’X)2'*p  X’X 


TTdS: 


N-j+i 


dX’X 


integrating  the  angles  out,  one  obtains  g(X*X)  dX’X,  mentioned 

above,  so  that  the  conditional  d ? str i but  ion  of  X  given  X’X  is 

t  (X  X’X) dX  =  7jd5:N_J+l 
J  ST 

m-j+S 

when  expressed  in  function  of  the  angular  coordinates  of  the  polar 
systems  of  coordinates. 

dS„ 

N-j+l 

is  the  distribution  of  a  unit  vector  Y-1  equ  id  i  str  i  but  ed  on  the  unit. 

sphere  of  the  N-j+l  dimensional  space. 

The  distribution  of  X  given  X’X  is  independent  from  the  nature 

of  f  ,  therefore,  al!  the  results  of  3 .4  and  3-5  are  valid  without 

the  assumption  of  normality,  it  is  enough  that  f (X)  =  h(X’X). 

making  the  change  of  variables  X  =  YT*  and  integrating  out  T’ ,  we 

obtain  the  distribution  above,  (since  T  and  X’X  are  equivalent  as 

sufficient  statistics),  i.e.  the  Y  have  the  distribution 

rTd^L  /5~  ;  more  specifically,  the  random  basis  (Y<  I  9  o  e  -y  YP) 

'j1  N-j+l  N-j+l 

is  constituted  with  unit  vectors  Yd  that  are  uniformly  distributed, 
given  Y  on  theunit  sphere  of  the  N-j+l  dimensional  space 

and  the  distribution  of  Y  is  the  same  as  that  of  X  given  X’X  =  J_e 
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Y  =  XT**"*  can  be  considered  as  a  student i zat Ion  or  normali¬ 
zation  of  X  generalizing  the  normalization  by  V  Nm*(X)  in  the 
univariate  situation.  One  can  write 

tt^-m  _  •nrc(N-j+i)/23dZM.1+1 
'j 1  sVj+i  '  j2PTrNp/2-p(p-')/4 


i  \*  j* 

YJ  has  coordinates  Y"  Y*:  with  respect  to  any  orthonormaf 

j  N 

basis  spanning  the  N-j+l  space  it  belongs  to,  and  the  distribution 
of  those  is 


r  [ (N-j+l )/2]dYJ 

J 


I# 


dY 


*J 
N- 1 


T 


(N-j+!)/2 


CMY 


j*2 


j*2  1/2 

N-  r  J 


2 1”*  t(N~j+S)/2)  dYj;c.0dYJ" 

'  "  j  N- i 

F{N-j+?)/2  — 7-7* 


(Y**  can  take  two  values  when  the  other  dependent  coordinates  of 

,N  - 

YJ  are  fixed.) 

3.6.2  To  obtain  the  distribution  of  X  given  X*X  as  an  explicit 
function  of  the  X^  or  the  Y^,  one  ought  to  pass  from  the  Y*  to 
the  Y,  choosing  Np-p(p+!)  functionally  independent  Y^,  then,  if 
needed,  come  back  to  X,  through  X  =  YT*j  only  NP-p(p+l)/2  X^ 

i 

are  functionally  independent,  given  X'X,  As  a  matter  of  fact, 
only  l-J-j  coordinates  of  X J  are  independent,  given  X*X;  this,  in 
view  of  the  fact  that  T  is  triangular  and  that  only  14— j  coordinates 
of  Y«i  are  independent.  This  shows  that  it  is  possible  to  get  the 


distribution  of  a  subsample  of  k  observat i ons  ? ,  given  X’X,  only 
if  the  minimum  of  N-j,  namely  N-p,  is  at  most  equal  to  k,  As  an 
explicit  function  of  the  independent  Y^  (that  we  will  denote  as  Y*} 

i 

the  probability  distribution  of  Y,  i.e.  of  Y*  is 

2PD!Viif±^  irdY’ 

D(r*)  TT^-M-j+i 

To  obtain  the  Jacobian,  one  can  refer  to  [10]  where  it  is  shown 
that  the  Jacobian  D(X)/D(T, Y":')  of  the  change  of  variables  X  =  YT? 
is  given  by 

P(X)  _  0(  X,  Y*  Y)  ,  D(  Y*  Y) 

D(T,Y»)  0(1, Y)  *  D ( Y*# ) 

(where  Y-^  denotes  t no  set  of  dependent  Yj*) 

Mow,  D(X,  Y{  Y)  __  p -rr*  ,  n~  t  see  [  10] 

D(T,Y)  ~  !  T'5 

Shd  ^  (  ~i~)  _  o"PTT 

D(X*  X)  (  if 

so  that 

(N-p- 1)/.? 

(X*  )dX  =  h ( X * X )  X;X  dX’X  D( Y»»)  dY* 

D(Y'Y) 

’  >  ^H-p+l  ^ 

D(Y*)  D(Y«Y) 

the  distribution  of  Yi;'  is 
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(2P/TTr  )  dY*/[D(Y’Y)/D(Y**] 

N-j+S 

N  |  t 

where  the  general  element  of  Y’Y  is  21  Y.YJ  and  of  Y*#  is 
£  5  "  h 

YN-k+l  *or  **  4,  *  an<*  °  *or  ^ 

3.6.3  In  case  N  =  p,  one  can  also  use  Cayley’s  parametric  repre¬ 
sents!  ion  of  Y. 

If  j 1+Y  |  #  o,  i.e.  if  Y  is  non-except i ona I ,  then  there  exists 
a  skew-symmetric  matrix  S  such  that 

Y  =  (i+sr'd-s)  =  (I-SHI+S)’8  =  (X-S)/(I+S) 

S=( I-Y) ( I+Y) ~  *  =  (I+Y)~S(r-Y)  «  (X-Y)/(X+Y), 


S  non-except iona I  since  skew  symmetric- 


if  jl+Yj  =  o,  i.e.  if  Y  is  exceptional,  then  it  is  known 
that  if  one  defines  a  sequence  ^ ^  of  diagonal  matrices 
(  +  !,.»„+  1} 

J ,  —  (IS  ....  I ) 

J2  =  (-!,  Il.ee.  i) 

J  =  (-1,-1,  l..u.  I  ) 

in  such  a  way  that  Jr  differs  from  by  changing  the  sign  of 

only  one  diaaonal  element  and  if  M.„  is  that  set  of  matrices  Y  such 
that 

Y+X  =  o 

|J2Y+I|  -  ° 


JJr-jY+I  j  =  ° 


J,.Y  -i-  X 


the  set  of  exceptional  matrices  is 


0 


Mr  and  Cayley’s  represen- 


i  as- 
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tat  ion  for  such  a  matrix  is 
Y  =  jr  (r-S)/(I+S),  S  =  (I-JrY)/(I+JrY) 
for  some  r>  I 

Let  i.i s,  then,  partition  the  space  of  X  into  the  sets  X*  and 
X"'*  such  that-  respect i vely, 

X*  ~  Y*T  Y*  n  o  n  -  ex  c  ep  f  i  o  n  a  5 
and  X*"  =  y#*T,**>  Y**  exceptions 3, 
one  has 

f  (X*)dX*  =  h(X-;:'*X*)dX* 
let  us  make  the  change  of  variables 

X*  =  [(  l-S*)/{  H-S“')  ]T{  ® 
whose  Jacobian,  after  Hsu  Eli]  is 

D(X*)/0(T*S*)  =lf t?7t  2?(p”!)y2Jl+S«|  ^ 


? 


i  I 


i  -i/2 


then 

dX*  =  2_P2P '  p"  ^  / 2j  X* » X*| 
hence,  the  distribution  of  S®  given  X**X*  is  independent  from 


~( p- ! ) 

i+S*  dX*fX*ds* 


X*fX*  end  given  by 

I  I  p~ 4  / 

KdS*/  J  I->-S‘"j  with  K  =  !  dS*7|I+S* 

'"3 

K  is  given  in 


p-i 


as 


,.(P_2)(P-0/21^irK/2/p(K(,2) 


I  n 


Co  c 


exceptional,  one  makes  the  transformat  ion 


Xs*  ~  J  t(I-S**)/r+S**)]T»**  when  Y**  £  M 

r  r 


It  results  from  [12]  that  <’A?  has  same  measure  as  M j  and  Mr  has 

measure  zero  if  r>  2.  Considering  only  the  set  we  will  have 
,  f  — i /2  i  |  — (  p—  ! ) 

dx**  =  2"p2P'P",i/‘-  X'**X»*j  j  I+S‘;::: j  dX •  **x**dS** 

and  the  distribution  of  S**  given  is  the  same  as  that  of 

S*  given  X»*X*. 


i 

I 

i 


s- 

I 

} 

ti 


-  !  7 


Given  the  proper  definition  of  on  M^,  S=S41"»‘'  on  we 

will  have 

f(X)dX  =  2“p2p(p‘S)/2hCX'X)(X»X)“S/2jl-t-s|  ^“^dX^XdS 


f .  e. 


and  given  X'X,  S  has  distribution 

[{2P(p-l)/2)/(^p  )]  dS/|l-i-s|  P“  ' 

£  P-J+i  1 

-p(  p+ 1 )  /4j^*  T  i  (  p-  j+ !  )  /2]2 ! (  P~  ' ]  ( P~2  }  t2~  ‘  |  I+sj 


-( p- 1 ) 


dS 


3*7 


Consider  +  he  space  defined  by  all  X  given  X*X,  if  G  is  a  group 
operating  transitively  on  xjx»x,  there  exists  at  most  one  probabil¬ 
ity  measure  on  that  space  that  is  invariant  under  G. 

if  we  consider  the  group  G  of  orthogonal  matrices^O^  operating 
on  X„  X-OX,  ft  is  geometrically  obvious  that  G  operates  transitively 
on  xJx'X,  on  the  other  hand,  it  is  straight  forward  that  the  condi¬ 
tional  distribution  of  X  given  X’X  is  invariant  under  G  since  (OX  )*(0X) 
X’X,  therefore  it  is  the  unique  distribution  on  xJx*X  invariant  under  G. 

Suppose  now  we  want  to  "pick  at  random"  a  p  dimensional  orthonormal 
basis  in  the  N  dimensional  space;  by  this  is  meant  choosing  X  such 
that  X*X  =  I  and  that  the  probability  distribution  of  X  be  invariant 
under  the  orthogonal  group;  it  sufficies  to  take  X  with  probability 
d istribut ion  f (X)dX  =  g(X*X)dX  and  normalize  by  X-YT' ,  Y  wf I  I  be  the 
basis  desired.  To  construct  such  a  basis,  one  can  take  X*,„.„,XP 
and  use  the  Schmidts*  orthogona I izat ion  process.  Alternatively,  one 
sees  that  one  can  choose  at  random  a  unit  vector  Y*  in  the  N  dimen- 

o 

sional  space,  then  a  random  unft  fector  Y  in  the  space  complement 
to  Y* ,  and  so  on;  at  the  jth  step,  one  takes  a  random  unit  vector 
In  the  space  complement  to  that  spanned  by  Y*...YJ"''*0  3y  random 
unit  vector  is  meant  a  vector  that  is  equidistributed  on  the  unit 
sphere  of  the  N-j+l  dimensional  space. 
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3.8  From  3.5  we  see  that,  in  case  f (X)dX  =  h(X»X)dX,  a  minimum 

variance  unbiased  estimate  of  a  function  of  the  parameters,  if 
it  exists  and  except  for  a  condition  of  completeness  on  S,  will  have 
to  be  found  among  the  solutions  of  the  integral  equation 

VS)=[---  Jfk(S)f  (S5  i  S,dS.  .. 

4.0  Further  generalizations  of  Thompsons  disfr ibutions  can  be 

tried  in  several  directions. 

4.!  One  can  try  to  obtain  the  conditional  distribution  of  a  sub¬ 

sample  5  of  a  sample  X  with  probability  distribution  -f  (X)dX, 
given  the  mean  Y  and  the  variance  S3  of  the  sample, in  the  general 
case, when  f (X)  is  not  spherical.  This  is  possible  when  one  has 
the  probability  distribution  of  the  couple  of  statistics  Y,Sa  for 
a  sample  of  any  size;  references  [5l>  [6],  [7]  deal  with  the 
derivation  of  such  a  d isfr ibut ion* 

l  * 

I  ••• 

and  of  the  variance  and  mean  S2,Xn  of  an  independent  sample  of 
size  n  -  N-k,  (where  M  is  the  size  of  X),  makes  the  change  of 
variables,  (as  shown  in  [2]), 

7n  “  7n/(N~I0  -  ?  k/(N-k) 

Ss  *  S2N/(N-k)  ~  (§  -  X)3/ka(N-k)3  -(?  -  7)' {Z  -  7)/(N-k) 
n  • 

whose  Jacobian  is  (N/N-k)  and  divides  by  the  distribution  of  X, 

2 

s  # 

If  (5Tn,Sa)  is  the  probability  density  of  one  has 


One  writes  the  joint  distribution  of 


5  -  (?■ 


f  (  £  |  7,S*)  =  tn^V7,  *  >!  Sn(5r’S“’  £  >Si  1](H/N-k)a>(  S'  )d? 

tN(5f,S») 

over  the  proper  domain  (namely  S®  >  o)  and  the  recurrence  formula 

Yuw, sa)  =  J.o  /^n(^  r >, s®(7ssa,  r  .<  >f (n{N/N-k)adi 

where  the  summation  is  performed  over  the  proper  domain. 

This  formula  may  be  useful  when  y  can  be  obtained  directly 
without  too  much  difficulty  for  sma S I  n. 

In  case  n=2  one  will  perform  the  rotation 

Yt  =  ( Xx-Xs)//2  Xx  =  ( Y  t+Y  a  )  /\fz 

Y2  *  (Xx+Xa)//2  =  1/27  X3  =  (Ya-Yi)//? 

7  =  Y 3//zT  S2  *  Y?/2o 

One  can  also  be  interested  in  the  distribution  of  the 
studentized  ¥  ,  namely  f  =  (5  -50/S, 

One  writes  the  joint  distribution  of  ?  =  (  £  xo „  and 

of  S2,>Tn,  makes  the  change  of  variables  t  =  (  £  -  ^n)/sn  and 
integrates  out  X^S^  the  density  of  t  is 

"<0  -  jf)WSn>J  <*„  +  «n>VVSn 

then  one  makes  the  change  of  variables 
t  =  (t  +  Tk/N-k) S/Sn 

333  I /2 

=  ( t+Tk/N~k)/[N/(N~k)~t  k  /(N-k)  -t*t/(N-k)] 

=  [t  +  Tk/(N-k)]/{(N/N-k)[l  -  T“k/(N-k)  -  S^k/N]$  UZ 


The  Jacobian  of  the  transformation  is 


=  C(N-k)/M](k_2)/2  \\  -  [T2k2/CH*k)  +  t* t]/N  ^  ”(k+2)/2 

«  [(N~k)/N](k~2)/2  0  -  T*k/(N-k)  -  S®k/wr(k+2)/2 

This  can  be  shown  as  follows,  let 

/=  CW-t *t  -  k273/(N-k)]/(N-k) 
y  =  t  +  k7/( N-k) 

direct  computations  shows  that 

gf ^  *3K/2(M-k)"kJc(M-k)^fj.  +  ni  +  if?  tf  j| 

the  determinant  is  that  of  the  matrix 

(N-k)  d  (I  +  [!/(N-k)]u*u  +  C  !/(N~k)  /  ]  if  *lf  ) 


where  u  is  a  row  of  one  u  =(!!,„ .1 1 ; ,  it  results  from  the  identity 


I  +  XxYx  +  X3Y2 


that  this  determinant  is 
k  n  k  ' 


1  +  YxXa'  YgX,1 
YjXg  I  +  YaXs 


(N-k)  l  C(i  +  W  )(|  +  k  ) 


ka?2  3 

Trap? 


*  (N-k) 


k-2  ^k-3  s 


which  gives  the  result  above.  Then,  g(f)dt  is  easily  obtained  as 
g(t)  =  h(t)D( t) 

EuT 

In  case  k=l  this  reduces  to 


I 


l- 


t. 
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t  =  t[N(N~!)]l/2/( l-fs/N~l) 1/2 
=  CN(M-!)]l/2/(!-ts/N-!)3/2 

4.3  Another  line  of  generalization  arises  from  the  consideration 

of  distributions  admitting  a  sufficient  statistic  (scalar  or 
vector).  The  problem  is  to  find  the  conditional  parameter  free 
distribution  of  a  subsample  ?  of  size  k  given  the  sufficient 
statistic  obtained  with  a  sample  X  of  size  N.  Nothing  very 
general  can  be  said  about  this  prob ! em  wt thout  making  further 
assumpt i ons. 

in  case  the  distribution  of  an  individual  observation 
belongs  to  Koopman*s  family,  with  one  parameter)  namely 

•f  ( X  f )  dX  j  =  exPCg(0)+t(X?)h(S)  +  4>  (Xf )  ]dXf ,  h*  *  0 


a  change  of  scale  of  the  parameter  in  an  interval  where  h(0)  is 
monotone  gives 


f  (Xf)dXj 


g(s)  2 


t(X;)8+/(Xj) 

dXf 


Then,  it  is  well  known  that  the  character f sr ?c  function  of  the 

II 

sufficient  statistic  T(X)  =  t(X.)  is 

i 


\D  (u)  =  [g(Q)3N/Cg(G+?u)]N 
1  T(X) 

This  is  straight  forward  since 


f  tut(X)  D  f  t(X)CFu+0]W(X) 

Je  t  (X)dX  =  g(Q)  =  Je  /  dx  =  g(G)/Cg(0+?u)D 
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> 


In  any  case,  the  probability  distribution  of  T(X)  can  be  obtained 

by  recurrence,  using  the  classical  convolution  formuleSo  Under 

r  14] 

proper  conditions  of  continuity  and  boundedness,  and  with 
obvious  notations, 


■WW  ■ 


pkCW 


pn(WVMV'1Tl< 


To  obtain  the  distribution  of  a  subsample  5  of  size  k  given 

the  sufficient  statistic  T(X)  of  a  sample  X  of  size  N,  one  gets 

the  distribution  of  T(  5 j  )  , ,  0T(  ?  k)  given  T(X)  and  comes  back 

to  the  original  variable  ?  „  The  characteristic  function  of 

T(  ? J  ) , „ J( ? k)  is  Cg(C) 3k/  TT  g(©+iuj) ,  that  of  Tn„,k  = 

j 

T(X)  -  T(5  )  is  C^PC©)/^  (©+iv)]N~k  that  of  T(  ? f  } .  0 ,T(  £  k)/n„k 
M  N  k  ^ 

is  [g(o)3  /[if  (o+i  v)3  JT  g(0+i uj);  hence  the  character ist ?c 

j 


function  of  T(  ?  „T(  Sf  ,}T(X)  is 

2  K  / 


[g(e)]N/Cg(©+iv)]N'k]T  gC0+l(Uj+v)3 

j 

UjT(?r  )  +  0+U.  T(¥  )+vT(X) 

since  E[e  !  5  k  k  3 

(u  j+v)T( ^ 5  )+<,+(  uk+v)T( £  k)+vTN_k 
=  E[e  3 


**  y-  .  -  •  '■'••'  '  *  •  *  V  -  -  -i  - 
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from  which  one  gets  the  characterist tc  function  of  T(  £  j) . 0T( £ 
given  T(X)  whose  charact er i st i c  function  is 

CgCS)/g(8+iv)]N 


The  result  is 

ip  fu,,  — ,* k) 

T(  ?  t)..T(  %k)  |  T(X) 

It  is  always  possible  to  obtain  the  probability  distribution 
of  ?T  given  T(X)  by  writing  the  joint  probability  distribution 
of  T(  ^...K  Ck)  and  T  N_k  =  T(X)  -  T(  S  ),  making  the  change 
of  variables  T^_k  =  T(X)  -  T(?  )  and  dividing  by  the  distribution  Qf 
T(X).  To  obtain  the  distribution  of  T(?  )  given  T(X)  one  writes 
the  joint  distribution  of  T(  ¥  )  and  T^,  makes  the  change  of 
variables  TfJ_k  =  T(X)  -  T(?  )  and  divides  by  the  distribution  of 
T(X),  If  T  =  T(X)  is  complete  and  if  f  ip  n(Tn)^  Is  a 
sequence  of  unbiased  estimates  of  a  function  g(G)  of  the  parameter 
Q,  then 

^(T)  =  J^kCT(?  )]f(T(t  }  |  T)dT(C  ) 


(e“lvT(X)[g{8+fV)3k“’NTT(g[G+i(u-+v)]rldv 

i _ _ j  _ J 

Je-fvT(X)j.g^s+fv-j-.N 


i.e.,  minimum  variance  unbiased  estimates  are  to  be  found  among 
the  solutions  of  the  integral  equation  above.  Also, 

{...  C  f  |  T)d  € 

will  be  an  unbiased  minimum  variance  estimate  of  E[\p(£)3 
if  if  exists. 


t 
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Since  the  character fst ic  function  XA(X|)  of  a  set  A  is  an 
unbiased  estimate  of 

the  best  estimate  of  this  measure  will  be  given  by 

and  consequently  the  best  estimate  of  §  (Xj)  by  ^ .  j  T  ) 0  Tne 
discrete  case  can  be  treated  in  a  similar  way. 

4«4  in  this  paragraph,  the  straightforward  technique  mentioned 

above  is  applied  to  several  important  d i str i but ions0 
a]  Poisson’s  distribution. 

In  case 

PCX- )  =  e“mmX‘-/X  » 

*  f 


X ,+.  o  +X..  =  N7  is  a  sufficient  statistic,  Poisson  distributed 
J  N 

and  it  is  complete,  one  obtains 

nT  kf  (NT)  i 

PC?  T)  =  [i-(k/N)]  [i/(N-k)3  P  „  -  k  *  N 

»  I*...  ?kJ[NT  -  kf]J 

on  the  domain  NT  -  k  5  >  o 

In  case  k=N,  one  obtains  the  well  known  result 

■  k  5  j 

p(?  MO  =  <N7)s/(TT  n  1  %  *  ) 

a  J 

the  best  estimate  of  P(X.  L  X)  is 


P(^f4x)  =  ^  ri  -  (I/M)]  L  I  /  ( N- 1 )  J  *(  6J)  t  X  4 


P(f|S)*(k/N>  [I  -  (k/N)](kj).  (N)T-k? 

This  is  a  Binomial  distribution, 
b)  Gamma  distribution  P( a,  h) 

f  (Xt)dX.  =  Cah/  P(h)3x^'  e’^'dXi, 

if  h  is  known  Y  is  sufficient  for  a,  complete,  and  N3T  is 
P (a,Nh)  distributed,  one  obtains 

f(?|7)dl  =  - - TT  (  ?  {/NX)h~' 

[P(h)]kr[(N-k)h]  1 

[!-  (  S’  /5T)(k/H)  3<H“k)h-l  d|  /(N3<)k  i  k  A  NY 


i.e.  the  £j/n7  are  parameter  free 
the  d  i  str  i  but  ion  of  £  j  5T  is 

h  !  i  k  r  kh-i  k  r  (N-k)h-l 

|(  S  *)d  ?  -  - - -(S  ^-)  (I-  —  — ) 

1  B[kh,  (N-k)h]  N  Y.  N  5T 


d(k  1/NY) 


i-e..  k  ?/M^  is  8[kh,(H-k)h]  distributed. 

The  best  estimate  of  P(Xj  A  x)  is,  for  x  A  N)?, 


^i.h-S  r.  (N-J)h-I 

■;-)  (i-  ii)  dX./M 


I,  A  NY 


b*)  in  case  h=l,  P (a,h)  is  the  exponential  distribution,  one 


obta ins, 


f  (  X  7)dl  =  C r*  CN)  /  P(M-k)]  .( I- S  il)N"k*’1  .d  ?/(N7) 


f(S  K)d  £  »  C 1/  3(k,N-!<)]  (-  2.)  (1-ai)  d(kf/NX) 

mV  mV 


k  |  £ 


The  best  estimate  of  P(X.  £.  x)  is,  for  n  ^  NT, 


Cl/  B(  i  ,N- 1 )  ]  (  !-  5{/N*.  )  d(  ?  ./NX)  £  /  N* 


b" )  The  case  of  the  sc-cailed  WeibuH’s  distribution 


/,  .  cc-i  ~{px.)a 

^(Xj)dXj  =  *P(PX{)  e  dXj 


for  known  a  is  not  different  from  that  of  the  exponential 
distribution  except  for  a  change  of  variable. 

Binomial  distribution 

Xj  is  0  with  probability  i-p  and  1  with  probability  p, 

X*  =  Xj  *.  „  is  sufficient,  complete  and  Binomial  disfri- 


but  ed 


pc  s  I  x»)  -  c”;^  -5  .)/<"„) 

with  =  £  1  +  ,,.  +  '5f; 


pi  !  x-')  =  (k  )(N_k  )/(k  ) 

k  *  '  1  x*-  x*; 

a  hypergeometr i c  disfribut ion. 


Multinomial  distribution 


X*  Xm 


P(X?.,.Xm)  =  (Nl/X*:...X  Up  1  ...pm  ,  X»+,00+Xm  a  N 

(X* , „ . . ,Xm)  is  a  complete  sufficient  statistic  X.  *.(  ?•..? 
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p(  ?|x)  =  (Jl)..(xm)/(J)  p£  4  xl 


V'**'  'k'  5 

e)  Rectangular  distribution, 

f  (X  . ) dX .  =  dXj/P,  o  <0  Xj  /.  9 j 9  >  o 
The  observations  being  ordered  X^^  ^  4:  * 

is  a  sufficient  and  complete  statistic. 


One  picks  a  set  cf  k  observations  ? ^ j x . . .  ? 


(k)* 


The  probability  that 


( k) 


Xf 


00  1 


is  k/IJ  and  that  X ,  k)  *  X(N)  is  [l-CK/M)], 


*(k)  #X(N) 


(X  X  ) 

{  s  (I)-*  Ok)1 


X(N))d^  =  k!C(H-k)/K]d?  /X!<N) 


fa 

i  (  &  (k)|  X(M))d  ^(k) 


°  ^  k )Z-  XN 


=  C(k(n-i<))/M3(?(k)/x(M))u-  dft(k)/x(N) 


k-l 


°  ^?(k)^XN 


•?}  Tps  t  ei  n-Sobel '  s  distribution 
p  -a(X--G) 

|  (X.)dX,  =  ae  ! 


On  i 


Xj  ^  G 


A  /  J  x  1*4  i  A  ^  ;  x  4  4  0  Y—  A  ^  |^j  ^ 


then  (Xfi'-,  X  -  Y.i,\)  is  a  sufficient  and  complete  statistic. 

The  probability  tnat  ?  =  X,  .  is  k/N  and  that  X  .  lStX.  . 

\  i )  UJ  t  I  /  \  I ) 

is  I  -  (  k/’M  )  . 

The  dist;  i but  ion  of  ( X ^  | ^  ,  X  -  =  Y)  is 


(aN)_ 


r  (m 


N-I) 


'N“2e  "'la^x(D  +  Y  ~  s)dX^dY 


cons i der i ng 
of  X  ,  then 


/-X.,.,  lei  X’  be  the  complementary  sample 
(  I)  (I) 
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k  * 


X(!)  =  X(0 

Y  =  7  -  X  =  (N~k)7*  +k  ?  -  X(}v  =  C(N-k)/N]Y»+(k/N)( * 

(I)  m  .  ' 

l.e.  Y*  =  [ NY/( N-k ) ]  ~  (k/N)(  ?  -X.) 


the  usual  technique  gives 


d  ?/  Nk(^-X(  ,j)k 

in  the  domain  where  the  quantity  between  brackets  is  positive* 

f(  $(!y  ?  “  ^(l)  j  x(j)»7  "  X(|))  = 

m  p  k(  f  ~  t\)  k“2  k(  C  -  X(  ,%)  N-k-2 

k(k~l)(N"2)[(N-k>/N:[ - Cl  -  - ^-3 

N(7-X(j))  N(7  -  X(I)) 

dk  £ (  | )  ^  dk(  £  -  ^ |  ) 

M(*-X(,))  N(7-X(1)) 

in  the  same  domain* 

The  first  formula  is  In  agreement  with  a  result  by  R«  F*  Tate, 
CI53,  P»  361,  formula  (?.8). 
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Log-Normal  D  IstrlbutFon 
Estimation  Problems 

by  Andre  G,  Laurent 
Department  of  Mathematics 
"/ayne  State  University 

The  use  of  the  Log-Normal  distribution  has  been  rarely  advo¬ 
cated  In  the  literature  devoted  to  industrial  life  testing  and 
failure  theory,  though  considerable  emphasis  has  been  put  on 
that  fool  In  papers  dealing  with  sensitivity  data  and  it  has 
been  widely  used  as  a  model  in  the  field  of  biological  assays, 

(Refer,  however,  to  "The  Exponential  Distribution  and  *ts  Role 
in  Life  Testing",  Technical  Report  No.  2,  May,  1958,  by  B.  Epstein). 

A  somewhat  extensive  exposition  of  the  history  of  the  Log- 
Normal  distribution  is  given  in  [I]. 

A  broad  variety  of  realistic  and  "reasonable"  hypotheses 
about  the  mechanism  of  failure  leads  to  a  Log-Normal  distribution 
for  the  age  of  death  or  the  intensity  of  the  stimulus  under  which 
failure  occurs.  Among  them,  that  class  of  hypotheses  which  involve 
random  effects  that  are  multiplicative  instead  of  additive  and 
i ead  to  the  use  of  the  centra!  limit  theorem.  For  example,  if 
failure  is  caused  by  an  accumulation  of  infinitesimal  random 
shocks  whose  effects  are  not  independent  but  are  sequentially 
proportional  to  the  already  accumulated  total  effect.  This  is  a 
well  known  set  up  and  a  plausible  model  for  failure  data  when  an 
ageing  process  takes  place. 


Let  X  be  Log-normal  a  I  atr  1  buted ,  i.e.  Y 

distributed.  The  moments  of  X  ere 

r2<j2/2  rm 
Mr(X)  -  e  e 

whence  the  expected  value  and  variance  of  X 


Log  X  is 


cr2/  2  m 


*  •  ^ 
f»l^  ~  U  U  ✓>  J  ~ 


ua  “  var[X]  =•  e2m(ec  -1) 


wh  ?  1  e 


E[Y3  =  m, logarf thmic  mean 


var[Yl  <=  o~.  logarithmic  variance 
i  i  n  i  rnun 
Co/(n-l)]s" 


rhc  minimum  variance  unbiased  estimates  ot  in  ur.d  c- '  being  Y  and 
.2 
V 


m  and  cr"  are  measures  ot  central  tendency  and  dispersion 
in  the  lognri thmic  scale,  useful  correspond Ing  characteristics 


in  the  original  scale  of  X  are  :i  =  em  that  is  the  median  of  the 


e  or 


distribution  and  ^  - 

"Naive"  estimates  cf  p  on. 


t  , 


I  2  <sy 

yi  a  re  e  ‘  and  e  , 


3 

Sv 


!/ n 


e  -  exp['g!  vog  X » ]  -  (TTX;)  is  the  geometric  mean  of  the  X  * 


These  two  estimates  are  not  unbiased; 

_  y  m  2 

since  Y  is  li( ni,o  /  n)  Hi.  e  j  -  em.t°  /2u) 

?  2 

and  since  the?  generating  function  of  Sc  is  (  !-2o  *f/n) 


-( n- I ) /2 


E[  ]  =  (  I  -  2o  ‘7ii) 


-( n- 1 )  fZ 
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We  want  to  obi  at n  estimates  without  bias 

Ij  e-w-  /*■"  _  t  r  t  <  \,N  /i.  •  ir— ‘-1'  \  k  ■ 

also 

rtk+(n.|)/23  2k 


.-0*2/2h  _  ^  j-  (,5)k/|<>;3[<r2k/(2n)K] 


ECS2k]  -  2k  — - t - 

Y  rC(n-!)/2] 


(<r  /ft  ) 


/k  v  rUn-D/23  _2k 

..  ,'Z  J  . .  —  jy  1 1>  an  unoiased  estimate  of  cr 

p[k+(n-!)/2] 


i herel ore 

P  f(n-0/2]  I  C  (-Dk  ml 


s2k 

Sy 


P  [k+(n~l)/2]22k 


i®  an 


unbiased  estimate  of  e”0""  ^2n  but  this  Is 

f -  o'  ,'p  r,>  _ 

12  '  /  >j L  '  jf  L(n-l )/2J  T 


(n^3)'v  *>Y/ 


v/here  J  denotes  the  Bessel  function  of  first  kind  of  order 
(n~3)/2. 

Since  Y  and  SY  ore  indeicndent  vor tables 

e7[2(n"3)/2/S^n'3)/2]  ^(0-0/2]  J(n_3)/2(Sy) 

is  an  unbiased  estimoto  of  em  =  since  Y  and  SY  are 

•  1 

sufficient  stat I st lcs; r t  is  a  minimum  variance  unbiased 
estimate  of  u. 

i 

2)  One  has 

K_  2  _  k  2k 

ehcr  =  X.  [h  /k!]c r 
K 

therefore 


P  C(n-!)/2jZ  Chk/k!jCnk/2k] 

k 

her  2 


is  an  un- 


2k 


biased  estimate  of  e 


now  this  is 
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r  [tn-O/S]  ...  2("-3)/a  Iu  D)/.tSy|/Shr;) 

(Sy/i^)ln-3)/:’- 

where  I  denotes  the  modified  Bessel  function  of  first  kind 
of  order  (n~3)/2. 

Since  Sv  is  a  sufficient  statistic,  the  minimum  variance 

1  ?  ? 

CP  r“  jj*  yp 

unbiesed  estimate  of  e  and  o  '  are  respectively 

rt(n-l)/2]  2(r“3)/g  Itr..o)/,o(3Y1/^)  ond 

(sy./^){n'3)/;’ 

r  r(n-D/2j  2(n"3)/2 _  J(n-3)/2(SY/") 

(SY  /n)(n'3)/2 

As  a  by-product,  one  v/i  I  !  obtain  the  minimum  variance  un¬ 
bio  sed  estimate  of  any  moment  Mr(><)  of  X. 

S  i  nee 

eC«rT]  * 

on  —  h(r  2 

depende  only/  Y,  ond  we  have  on  unbiased  ect  irnate  of  e  that 

depends  only  on  3V,  making  h  -  (n-l)r  /2n  ,we  get  o  minimum 
variance  unbiased  astir.-,  to  of  >\r  az 

*rY.PC(n-l)/?.J_  i(n-3)/?.(rSY/irf  ) 

'  rr,(n-3>« 

Ref  erenoes 

[!j  J.  C.  Aitchinson  anci  J.  A.  Brown,  "The  Los-Normal  Distri¬ 


bution  , 


Combr idgs, 


WAYNE  STATE  UNIVERSITY 
Technical  Report  Distribution  List 
Contract  Nonr-2575(00) 


Naval  Activities 

Head,  Statistics  Branch  2 

Office  of  Naval  Research 
Washington  25,  D.C. 

Commanding  Office  2 

Office  of  Naval  Research  Branch 
Office 
3^6  Broadway 
New  York  13,  New  York 

Commanding  Officer  2 

Office  of  Naval  Research  Branch 
Office 
Navy  Noo  100 
Fleet  Tost  Office 
New  York,  New  York 

Commanding  Officer  I 

Office  of  Naval  Research 
1000  Geary  Street 
San  Francisco  9,  California 

Technical  Information  Center  6 
Naval  Research  Laboratory 
Washington  25,  D.C. 

Chief  of  Naval  Materiel  1 

Code  M533,  Room  2236 
Main  Navy  Building 
Department  of  rhe  Navy 
Washington  25,  D.C. 

Dr.  Craig  A,  Magwire  I 

Department  of  Math,  and  Mech. 

U,  S.  Naval  Postgraduate  School 
Monterey,  California 

Dr.  Mitchell  L.  Cotton  I 

Department  of  Electronics 
U,  S,  Navai  Postgraduate  School 
Monterey,  California 

Chief,  Bureau  of  Ordnance  (Ad3)  I 
Department  of  the  Navy 
Washington  25,  D.C. 

Chief,  Bureau  of  Ordnance  i 

Department  of  the  Navy 
Washington  25,  D.C. 

Attn:  A.  Rothstein  Code  ReUg 


Chief,  Bureau  of  Aeronautics  I 
Department  of  the  Navy 
Washington  25>  D.C. 

Attn:  Wualftv  Control  Dfv. 

Commander 

U.S.  Navai  Air  Missile  Test  Ctr. 
Folnt  Mugu,  California 
Attn:  Chief  Scientist 

Commander  I 

U,  S.  Naval  Air  Development  Ctr. 
Johnsville,  Pennsylvania 

Mr.  H.  R.  Thcman  I 

Bureau  of  Aeronautics,  AV-4422 
Department  of  the  Navy 
Room  ^796,  w  Building 
Washington  25,  D.C. 

Mr .  Francis  A.  Thompson  I 

Bureau  of  Aeronautics,  AV-3102 
Department  of  the  NaVy 
Room  IW63,  Building 
Washington  25,  D.C. 

Chief,  Bureau  of  Yards  and  Docks 
Department  of  the  Navy 
Washington  25,  D.C. 

Attn:  Mr.  ’Vo  I  man,  Code 

Quality  Control  Division  (ycc)  1 
Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  r'.C. 

Statistics  Branch  I 

Qc  Division  Qc5 
Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.C. 

Mr.  R.  E.  Wiley  I 

Mater I  els  Branch 
Office  of  Nava!  Research 
Washington  25,  D.C. 

Air,  Masao  Yoshitsu  ! 

Naval  Inspector  of  Ordnance 
Naval  Gun  Factory 
Washington  25,  D.C. 


2 


Miss  Besse  B.  Day  I 

Bureau  of  Jhips,  Code  3I0B 
Department  of  the  Navy 
Washington  25,  D.C. 

Mr.  A.  Lieberrnan  I 

Bureau  of  Ships,  Code  3730 
Department  of  the  Navy 
Washington  25,  D.C. 

Mr.  A,  S.  Mart hens  I 

Bureau  of  Ships,  Code  373A 
Department  of  the  Navy 
Washington  25,  D.C. 

Mr.  F.  Brown  I 

Bureau  of  Ships,  Code  363B 
Department  of  the  Navy 
Washington  25,  D.C, 

Mr.  H.  Weingarten  I 

Bureau  of  Ships,  Code  280 
D-->rtment  of  the  Navy 
VL  ington  25,  D.  C. 

Mr.  F.  R.  Del  Prior e  I 

U.  S.  Naval  Engineering  Experiment 
Station 

Annapolis,  Maryland 

Captain  6„  L.  Lubelsky,  USN  I 

Quality  Eva 'nation  Laboratory 
U«  S.  Nava!  Ammunition  Depot 
Crane,  Indiana 

Commanding  Officer  I 

U.  S.  Navy  /nine  Defense  Laboratory 
Panama  City,  Florida 
Attn:  Mr.  J.  Boyd 

Commander  I 

Materiel  Laboratory 
New  York  Naval  Shipyard 
Naval  Base 
Brooklyn  !,  Mew  York 
Attn:  A.  Walner 

Commanding  Officer  I 

U.  S.  Nava!  Radiological  Defense 
Laboratory 

San  Francisco,  California 
Attn:  Miss  M.  Sandomire 

Dr.  Jul fus  Lieblein  I 

Applied  Mathematics  Lab.  Code  820 
David  Taylor  Model  Basin 
Washington  7,  D.C, 


Air.  E.  J.  Nucci  I 

Bureau  of  Ships,  Code  819 
Department  of  the  Navy 
Washington  25,  D.C. 

Office,  Secretary  of  Defense 

Inspection  and  Quality  Control  i 
Division 

Office,  Asst.  Secretary  of  Defense 
(S  and  L) 

Washington  25,  D.C. 

Attn:  Mr.  Irving  B.  Altman 

Inspection  and  Quality  Control  ! 
Division 

Office,  Asst.  Secretary  of  Defense 
(S  and  L) 

Washington  25,  D.C. 

Attn:  A\r.  John  J.  Rfordan 

Office,  Asst.  Secretary  of  i 

Defense  (R  and  E) 

Room  3^984,  The  Fentagon 
Washington  25,  D.C. 

Attn:  Mr.  R.  H.  DeWitt 

Office  of  Guided  Missiles  I 

Office,  Asst.  Secretary  of  Defense 
(R  and  E) 

Washington  25,  D.C. 

Attn:  Mr.  Carlton  M.  Beyer 

Technical  Library  I 

Office,  Asst.  Secretary  of  Defense 
(R  and  E) 

Room  3EIO65,  the  Pentagon 
Washington  25,  D.  C. 

Army  Act i vi t i es 


Mr.  Si  las  W? II  jams,  Jr.  I 

Standards  Branch 
Procurement  Division 
DCS  Logistics,  U.S.  Army 
Washington  25,  D.C. 

Dr.  L.  S.  Gephart  • 

Office  of  Ordnance  Research 
Box  CM,  Duke  Station 
Durham,  North  Carolina 

Air  Force  Activities 

LTCOL  W.  C.  Marcuis  I 

Air  Research  Development  Command 
Baltimore,  Maryland 
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Mr.  R.  Riedenbender  I 

Hq„,  Air  Materiel  Command  (MCQFF) 
'vr Ight-Fatterson  Air  Force  Base 
Ohio 

Dr.  J.  A.  Greenwood 

Hq.,  U.  S,  Air  Force,  AFC1N-3BX2 

'Vashington  25j  D.C. 

M I scellaneous  Government  Agencies 


ASHA  Documents  Service  Center  5 
Knott  Bldg. 

Dayton  2,  Ohio 


Office  of  Technical  Services  I 

Department  of  Commerce 
'Vashington  25,  D.C. 

A\r.  Leon  Gilford,  ADSS  I 

U.  S.  Census  Bureau 
Washington  25,  D.C* 

Dr.  Joseph  Daly  I 

U.  S.  Census  Bureau 
'Vashington  25,  D.C. 

Dr.  Marvin  Zelen  I 


Statistical  Engn’g  Laboratory 
National  Bureau  of  Standards 
Washington  25»  D.C, 

Business  Concerns 

Mr.  I.  Be  I  son  I 

Mine  Safety  Appliances  Co. 

Ga I  I ery,  Pennsy I  van i a 

Mr.  John  K.  Borkman  I 

Whippany  Laboratory 
Bell  Telephone  Laboratories 
Whippany,  New  Jersey 

Mr.  W.  H.  Ciatworfhy  I 

Bet  t i s  Plant 

”'est  i nghouse  Electric  Corp. 

Box  1468 

Pittsburgh  30,  Pennsylvania 

Mr.  H.  F.  Dodge  I 

Bell  Telephone  Laborator r es,  Inc. 
463  West  Street 
New  York  14,  New  York 

Mr.  S.  Weiner  I 

ARINC  Research  Corporation 
1700  K  Street,  N„\V0 
Washington  6,  D.C. 


Dr.  G,  Kaskey  I 

Remington  Rand  UNIVAC  Div.  of 
Sperry  Rand  Corporation 
1900  'Vest  Allegheny  Avenue 
Fhi laddfchia,  Pennsylvania 

Dr.  John  F.  Hofmann  ! 

Systems  Laboratories  Corp. 

15016  Bentura  8lvd. 

Sherman  Oaks,  California 

Mr.  J.  C.  Duryea  _  I 

Senior  Components  Engineer 
Sperry  Gyroscope  Company 
Great  Neck,  New  York 

Un  ?  vers?  t i es 

Professor  Cyrus  Dorman  I 

Depart,  of  Industrial  Engn’g 
Columbia  University 
New  York  27,  New  York 

Professor  Acheson  J.  Duncan  I 

Dept,  of  Industrial  Engn’g 
The  Johns  Hopkins  University 
Baltimore  18,  AAaryland 

Frofessor  Eugene  L.  Grant  I 

Civil  Engineering  Department 
Stanford  University 
Stanford,  California 

Dr.  Eugene  Lukacs  ^  I 

Department  of  Mathematics 
Catholic  University 
’Vashington  17,  D.C, 

Frofessor  E.  G.  Olds  I 

Department  of  Mathematics 
Carnegie  Institute  of  Technology 
Pittsburgh  13,  Pennsylvania 

Dr.  Max  A.  Woodbury  I 

Department  of  Mathematics 
College  of  Engineering 
New  York  University 
New  York  53,  New  York 

Frofessor  Ralph  A.  Bradley  I 

Department  of  Mathemat i cs^ 
Virginia  Foiytechnic  Institute 
Blacksburg,  Virginia 

Frof  essor  El  I  Is  Ott  I 

Department  of  Mathematics 
Rutgers  Unrversi fy 
New  Brunswick,  New  Jersey 


Foreign  Addresses 


Professor  Maurice  H.  Beiz 
University  of  Melbourne 
Carlton  N.  3 
Viet  ria,  AUSTRALIA 


Professor  Tosio  Kitagawa 
Mathematical  Institute 
Faculty  of  Science 
Kyushu  University 
Fukuoka,  JAFAN 


1 


I 


Mr «  Cesar eo  Vi 1 1 egas  I 

Instituto  cfe  Asthenia  "tea  y  Estadistica 
Av.  J.  Herrera  y  Reiss  I g  565 
Montevideo,  URUGUAY 


Mr,  f.  F.  ’7a de,  statistician 
Ajuminum  Company  of  Canada,  Ltd. 
Kingston,  Ontario,  CA IADA 

Professor  0.  a.  S.  Fraser 
Department  of  Math  etna-  i  cs 
University  of  Toronto 
Toronto  5,  CANADA 


